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The 7+ jet production in p—p collisions at high energies is studied. Double-spin 
asymmetry All of the process is calculated by using Monte Carlo code SPHINX. 
The predictions for All of the 7 + jet production in the p — p collisions using the 
spin-dependent gluon distributions with positive and negative sign of AG(x,Q 2 ) 
are made at RHIC energies. A possibility to extract the spin-dependent gluon 
distribution AG(x) at STAR experiment and to discriminate between different 
theoretical scenarios is discussed. 
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1 Introduction 

The accelerator complex of polarized protons (RHIC-spin) at Brookhaven National 
Laboratory will be first proton-proton collider accelerating proton beam up to energy 
of 250 GeV with polarization ~ 70% and luminosity up to L ~ 2 ■ 10 32 cm~ 2 s~ 1 at 
y/s = 500 GeV. RHIC-spin will have unique capability to study polarization phenomena 
such as direct 7, jet, dijet, Drell-Yan lepton pairs and W ± , Z° Jl], Q, || production and to 
perform the test of perturbative and non-perturbative QCD. The main goal is to obtain 
direct information on the spin- dependent gluon and sea quark distribution and to clarify 
nature of proton spin. 

In our previous works ||, [|, we studied the double-spin asymmetry All of jet, dijet 
and direct-7 production in p — p collisions at RHIC energies as a function of sign and 
shape of spin-dependent gluon distribution AG and made predictions of All for the STAR 
experiment planed at RHIC || H . We would like to note that global fits J7| to the present 
inclusive DIS data cannot even fix the sign of AG, not mentioning its magnitude M. For 
example, one of the analysis favors a positive sign, while the recent NLO analysis [|Kj] 
yields negative first moment of the gluon density. Both positive and negative values of 
the sign of AG(x, Q 2 ) over a wide kinematic range (10~ 3 < x < 1) were considered in fTT 



The possibility to draw conclusions on the sign of the spin-dependent gluon distribution, 
AG(x, Q 2 ), from existing polarized DIS data have been studied in fl2f . The result of the 



DIS data analysis [12| on g™ supports the conclusion that the sign of AG(x, Q 2 ) should 



be positive. Also, the recent HERMES result [13] indicates that AG/G is positive in 
the intermediate x region. Nevertheless, the additional confirmations on sign of AG are 
required and the experiments for direct measuring of the AG/G are necessary. We also 
consider that the NLO QCD analysis with negative sign of AG(x, Q 2 ) of DIS data on 
g\ structure function could give additional constraints of quark and gluon distributions 
moreover usually used 0. 

The measurement of the 7+ jet production in polarized proton-proton collisions is more 
perspective for determining AG(x) because the reaction is dominated by the Compton 
Q9 ~^ 7? scattering and we can assume that Aq(x)/q(x) is known from polarized DIS 
experiments |T5|-[p^[. Some contributions to the cross section of 7 + jet production 



come from non-direct mechanisms such as qq annihilation, photon fragmentation and k± 
smearing Q. 

The main goal of this work is to made predictions of double-spin asymmetry All 
of the p + p — > 7 + jet + X process and to study the possibility of direct extraction 
of spin dependent gluon distribution AG(x) from this reaction at STAR detector at 
RHIC. The determination of AG by means of 7 + jet asymmetry is connected with large 
background, similar as in the case of prompt photon production. Two closely spaced 
photons resulting from decay of high pt neutral mesons (tt°, 77) can produce fake 'direct 
photon signal' in detector. Because the cross-sections of the main subprocesses for 'direct' 
mesons production in the pp collision (qq — > qq, qg — > qg, gg — » gg) are much higher than 
the cross-section for Compton scattering, we can expect the large yield of background 
(direct meson + jet) events. So, the additional aim of the present paper is to study the 

3 The discussion on the problem can be found in H and references therein. 



kinematic properties of signal (S) and background (B) and to find the corresponding 
algorithm for background reduction. 

The paper is organized as follows. The different sets of spin-dependent parton distribu- 
tions used in Monte Carlo (MC) simulations are described in Section 2. The procedure for 
calculation of All and extraction of AG(x), together with the results of MC simulations 
at -y/s = 200 GeV are given in Section 3. Conclusions are summarized in Section 4. 

2 Spin-dependent gluon distribution 



The 3 sets of spin-dependent parton distributions [[T1J and [TJJ are used to calculate 
the 7 + jet asymmetry All. First set is based on work of Altarelli and Stirling [|14[ 
and include a scenario with large gluon polarization AG. Second and third ones have 
been obtained by the phenomenological method |TT| including some constraints on the 
signs of valence and sea quark distributions, taking into account the axial gluon anomaly 
and utilizes results on integral quark contributions to the nucleon spin. Based on the 
analysis of experimental data on deep inelastic structure function gi the parameterization 
of spin-dependent parton distributions for both positive and negative sign of AG have 
been constructed. It was found that the Q 2 evolution of structure function gi is sensitive 
to sign of AG. We would like to note the both sets of distributions describe experimental 
data very reasonable. We shall denote AG >0 and AG <0 sets of spin-dependent parton 
distributions obtained in [|ll|] with positive and negative sign of AG, respectively. It was 
shown in |23| that phenomenological method reproduces the main features of the NLO 
QCD Q 2 -evolution of proton, deuteron and neutron structure function g\. Therefore 
the constructed spin-dependent quark and gluon distributions can be reasonably used to 
study the asymmetry A L l of 7 + jet production in p — p collisions too. 

Figure 1(a) shows the dependence of the ratio AG(x,Q 2 )/G(x,Q 2 ) on x at Q 2 = 
50 (GeV) 2 for gluon distributions |14| and AG >0 [[II] . The curve ]14| increases with x 



and then decreases to zero. The behavior is due to the structure of the ratio AG/G ~ 
x a (l — x)P . The monotonous increase of an second curve with x connects to other 
asymptotic of the ratio AG/G ~ x s at x — ► 1. Figure 1(b) shows the dependence of the 
ratio AG(x,Q 2 )/G(x,Q 2 ) on x at Q 2 = 50 (GeV) 2 for gluon distributions AG >0 and 
AG <0 fTJ. 

The integral contributions of gluons Ag(Q 2 ) = / AG(x,Q 2 )dx to the proton's spin 
are found to be Ag >0 ~ 2, Ag <0 ~ -3.4 at Q 2 = 10 (GeV) 2 JTJ, and Ag ~ 3 at Q 2 = 
10 (GeV) 2 HU. Because integral quark contributions to the proton's spin is practically 
the same for all the parameterizations, the difference in asymmetry All should reflects 
the difference in shape of gluon distribution (and sign of AG) used. We would like to 
note that large value of gluon polarization corresponds to relatively small contribution 
of qq —>■ gj process to the final 7 + jet asymmetry. Because this process has the same 
final state as a signal, it is the background for determining the AG. But, our calculations 
shows that, at RHIC energies, the cross-section for qq — > (77 process is practically for 
an order of magnitude smaller than the cross-section for qg — > jq process, for all spin- 
dependent PDF's used, so such a kind of background could be neglected, especially if we 
compare it with direct mesons yields. 



3 Monte Carlo simulation results 

For detailed study of the p+p — > j+jet+X process we have used the Monte Carlo code 
SPHINX |23| which is 'polarized' version of PYTHIA |25| . Jet reconstruction was done by 



the JETSET-subroutine LUCELL [25]. This routine defines jets in the two-dimensional 
(rj, 0)-plane, r\ being pseudorapidity and </> the azimuthal angle. STAR detector covers 
full space in azimuth and pseudorapidity region — 1 < n < 2. 

In order to have segmentation expected at STAR (Ar/xA0 = 0.1x0.1) we used 30 
77-bins and 60 0-bins in our calculation procedure. The values of LUCELL-subroutine 
parameters E± u and R = \/An 2 + A<p 2 were Ef a = 1.5 GeV and R = 0.7. 

The expected resolution of the STAR Electromagnetic Calorimeter AE/E~Q.l6/\/E 



was also taken into account [26]. To obtain total rates (Rate = a ■ L) of 7 + jet events 
at STAR we taken into account designed integrated luminosity at RHIC: L = 320 pb~ x 
at \fs = 200 GeV. The expected values of polarization of proton beams at RHIC 
(Pbi = Pb2 = P = 0.7) were also used in the analysis. 

All of the SPHINX calculations were done with a cut on the partonic transverse mo- 
mentum pt > 10 GeV jc. The simulations include initial-state radiation, final-state 
radiation, fragmentation and particle decay effects. 

3.1 Identification of 7 + jet events and background reduction 

After generating the signal event sample we employ the following cuts for reconstruc- 
tion of candidate 7 + jet event: 

• there is exactly 1 reconstructed jet with —0.3 < r]j et < 1.3, E? > 10 GeV 

• and there is a 1 photon satisfying — 1 < r/ 7 < 2, p? > 10 GeV jc. 

We would like to note that study of the background enter at this level of the analysis. 
The direct meson background was simulated by using all possible 2 — > 2 partonic hard 
scattering processes. The event is treated as a candidate direct photon event if there is a 
high px {jpx > 10 GeV jc) neutral meson (71*°, rj) in coincidence with the 1 reconstructed 
jet (E T e > 10 GeV) satisfying the pseudorapidity cuts defined above. 

Figure 2 shows the comparison of the direct meson + jet to direct photon + jet yields 
at STAR in p — p collisions at y/s = 200 GeV as a function of transverse momentum p? 
for the spin-dependent parton distribution AG >0 [|TTJ. The signal/background ratio is 
found to be S/B ~ 1/10. As seen from Figure 2 it is necessary to impose additional cuts 
for background suppression. 

In such cases, the so-called 'isolation cuts' method, developed by UA2 Collaboration 



27fl , have been used for meson background suppression. The candidate direct photon was 
treated as 'isolated' if no other photons or charged particles fall within a cone of radius 
Risoi around the photon. The idea is that direct meson usually has much more particles 
in its 'neighborhood' than direct photon. 

We tried to make more sophisticated algorithm for photon isolation comparing the 
energy and multiplicities of particles falling within the isolation cone for different values 



of Risoi- The results are shown in Figure 3. We have used 3 different values of isolation 
cone radius: Ri SO i = 0.26 (UA2 condition), R iso i = 0.5 and Ri SO i = 0.7. 

As seen from Figure 3 particles accompanying the direct photon are softer than those 
for the direct meson. Also, the use of larger cone radius allow better reduction of back- 
ground without significant loss of signal. So, the conclusion is that the combination of 
cuts on multiplicity and energy will give better background suppression. Our results show 
that highest gamma efficiency P 7 / 7 and meson rejection 1 — P 7 / 7r o (here Py/^o is meson 
contamination) values are obtained under following conditions: the isolation cone radius 
is Risoi = 0.7 and summed energy of particles (if there are any) around a direct photon 
candidate must be lower than 1 GeV . In such a case, 90% of mesons were rejected with 
gamma efficiency P 7 / 7 of 88%. 

Figure 4 shows comparison of the direct meson + jet to direct photon + jet yields at 
STAR in p — p collisions at y/s = 200 GeV as a function of transverse momentum pt for 
spin-dependent parton distribution AG >0 |TJ| with the isolation cuts imposed. 



Additional event-by-event discrimination between single photons and photon pairs 
(from neutral mesons decays) at STAR is based on measurement of transverse profile 
of the electromagnetic shower in Shower Maximum Detector (SMD). It is expected that 
this approach will rejects ~ 80% of photon pairs from meson decays with single gamma 



efficiency of ~ 80% |8L 29 



Figure 5 shows comparison of the direct meson + jet to direct photon + jet yields at 
STAR in p — p collisions at y/s = 200 GeV as a function of transverse momentum p T for 



spin-dependent parton distribution AG [ill], after the additional meson/ '7 discrimina 



tion with the help of Shower Maximum Detector. 

Finally, we can conclude that only combination of isolation cuts and meson/'j dis- 
crimination with the help of SMD will be able to reduce the background in the whole px 
region of interest for measuring of double-spin asymmetry All and extracting of AG(x) 
at STAR. 

3.2 Asymmetry of 7 + jet production 

In Monte Carlo simulations, the double spin asymmetry is defined through the differ- 
ence of cross-sections (numbers of 7 + jet events) for antiparallel (||) and parallel (ft) 
spins of colliding protons: 

A ° U -° n 1 N? +jet -N? +jet 

- U + - TT P b ANH Jet + NH Jet ' U 

with the statistical error: 
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where P&i and p>2 are values of the polarization of proton beams. In our case, only 
the 7 + jet events that passed all the above cuts were used for calculation of All and 
extraction of AG(x). 



Figure 6 shows the asymmetry An of 7 + jet production in polarized pp collisions 
at y/s = 200 GeV for 3 different sets of spin-dependent PDF's (AG >0 Q, AG <0 Q 



and |L4j) as a function of transverse momentum p?. The absolute value of An increases 
with px for all spin-dependent PDF's used. There is a clear difference in An between 
positive and negative sign of AG as well as between PDF's AG >0 [|l J] and [I 1] in the 
region p T = (10 — 25) GeV/c. 

3.3 Partonic kinematics reconstruction 

The Bjorken x- values of colliding partons is reconstructed by means of measured 7 
and jet transverse momenta p? and pseudorapidities: 

Xl = 7i — 2 — ' (3) 

2pT (e-'h + e -vj«) 

X2 = 7i 2 • (4) 



The usual procedure |2g, |31| is to associate the larger of the x\, X2 values with the quark 

x q = max[ii,i 2 ], (5) 

and the smaller one with the gluon momentum fraction 

x g = mm[xi,x 2 ], (6) 

with additional condition: 

max[xi,X2] > 0.2. (7) 

Comparison of the simulated to reconstructed values of the initial-state parton kine- 
matics for the events that passed all above cuts is shown in Figure 7. 

3.4 Extraction of AG(x) 

The main formula for extraction of AG(x) from 7 + jet coincidence events can be 



written in the following form |28 



A G(XgJ - r> r> ^^ 



Nl l - M T 



PuPm J2ti K IS ( X ^ Q*ya LL (cos9 t )/G(x gi , Q?) ' 

where N^ and N^ are number of 7 + jet events for reconstructed value of gluon mo- 
mentum fraction x g , for antiparallel (|l) and parallel (||) spins of colliding protons, 
respectively and N s = N^ + N^ . Here dz,£(cos0i) is the partonic spin-correlation pa- 
rameter for qg — > qj process, 9i being the scattering angle in the partonic center-of-mass 
(pCM) system. Assuming collinear partonic collisions, the absolute value of cos 9i can be 
obtained as: 

cos^j = tanh(±r]j), (9) 



where 

Vi = ±^~Vjet)- ( 10 ) 

The value of the partonic spin-correlation parameter an(cos6i) can be calculated in the 
LO pQCD |34|] and for qg — >■ qj process we have: 

"2 ~ 2 

a L L = , 2 . 2 , (11) 

s z + ir 

where s and u are the Mandelstam variables for the hard scattering: 

s = X1X2S, (12) 

u = -|(l + cos0i). (13) 

The factor Afp S (x q i, Qf) is a parameterization of the A\ structure function for the 
proton obtained by fitting the results of polarized DIS experiments. The value x q i is the 
quark momentum fraction and Q\ is momentum transfer scale, assumed to be Qf = p\j1- 
For the calculations done here, we used the following parameterization of A® p IS |32| : 

4fW Q 2 i)~< IS M = 0-01902 + x/ 01163 (l - e" 1 - 845 "*). (14) 

The parameterization of unpolarized gluon structure function G(x g i,Qf), evaluated at 
reconstructed value of x g and evolved to the value of Q^, was taken from [p3| . 

The results of extraction procedure for spin-dependent parton distribution AG >0 [[Ll[ 
at colliding energy of \fs = 200 GeV is shown in Figure 8. In order to discuss the quality 
of extraction procedure, we show in Figure 8 the comparison of the input polarized gluon 
distribution to values extracted from the simulated data as a function of the reconstructed 
Xgiuon- The input polarized gluon distribution is evolved to the most probable value 
Q2 = 50 {GeV/cf. 

In general, theoretical input and the values extracted from the data agree very well. 
The agreement with results obtained in || is also reasonable. The systematic discrepancy 
for Xgi uon < 0.08 reflects the ambiguity in assigning Xip) values to gluon (quark) in the 
small x region and ambiguity in the choice of the sign for cos 9i. Also, the behavior of 
extracted polarized gluon distribution at higher x-values reflects the role of condition 
(7) and results in relatively small change of the shape of the distribution for x g i uon > 
0.2. Nevertheless, this result shows that AG(x) can be extracted at STAR from the 
measurements of 7 + jet production in p — p collisions. At s/s = 200 GeV, STAR 
detector will give us the full information about AG(x) in the region x g i uon = (0.03 — 0.3). 

Figure 9 shows reconstructed values of xAG(x) versus reconstructed values of x g i uon 
for 3 different sets of spin- dependent PDF's (AG >0 Q, AG <0 Q and PR ). One can see 



that measurements at STAR will be able to discriminate between scenarios with positive 
and negative sign of AG. In the case of positive sign of AG, extraction of AG(x) at 
y/s = 200 GeV will give us possibility to distinguish between theoretical spin-dependent 
PDF's AG >0 [[IJ and [0 in the region x gluon = (0.03 - 0.2). 



4 Conclusions 

The 7 + jet production in p — p collisions at RHIC energies was studied. Double- 
spin asymmetry All of the process is calculated by using Monte Carlo code SPHINX 
taking into account parameters of the STAR detector. Special attention was paid on the 
procedure for suppression of meson background. The algorithm for isolation of direct 
gamma based on 'isolation cuts' method was developed. A procedure for extraction of 
AG(x) from simulated data sample was described in detail. The obtained results show 
that highest gamma efficiency and meson rejection values are obtained under following 
conditions: the isolation cone radius is Ri SO i = 0.7 and summed energy of particles 
(if there are any) around a direct photon candidate must be lower than 1 GeV. In 
such a case, 90% of mesons were rejected with gamma efficiency of 88%. Also, we 
concluded that only combination of isolation cuts and meson/7 discrimination with the 
help of Shower Maximum Detector will be able to reduce the background in the whole pt 
region of interest for measuring of double-spin asymmetry All and extracting of the spin- 
dependent gluon distribution AG(x) in the region x g i uon = (0.03, 0.3) at y^ = 200 GeV 
at STAR. 

Thus, the obtained results allow us to hope that direct experimental information on 
a sign of the spin- dependent gluon distribution will give constraints on the quark and 
gluon distributions and to clarify a role of quark and gluon orbital momenta in proton's 
spin composition p: 
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Figure 1 

The ratio AG(x, Q 2 )/G(x, Q 2 ) of polarized and unpolarized gluon distributions as a 

function of x for two different parameterizations (a) AG >0 P] , |[I3 ] and (b) AG >0 (11 

AG <0 H at Q 2 =50 GeV 2 . 
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Figure 2. 

Comparison of the direct meson + jet to direct photon + jet yields at STAR in p — p 
collisions at y/s = 200 GeV as a function of transverse momentum p? for spin-dependent 
parton distribution AG >0 [ ITf . 
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Figure 3. 

The multiplicity and summed energy distributions of particles falling within an iso- 
lation cone of radius R around a direct photon (left column) or a direct meson (right 
column). 
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Figure 4. 

Comparison of the direct meson + jet to direct photon + jet yields at STAR in p — p 
collisions at ^/s = 200 GeV as a function of transverse momentum p? for spin-dependent 
parton distribution AG >0 |IT| , with the isolation cuts imposed (see text). 
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Figure 5. 

Comparison of the direct meson + jet to direct photon + jet yields at STAR in p — p 
collisions at y/s = 200 GeV as a function of transverse momentum px for spin-dependent 
parton distribution AG >0 [J I], after the additional meson/7 discrimination with the help 



of Shower Maximum Detector. 



13 



-A A— _A— 



Vs - 200 GeV 
-1 < v < 2 
L = 320 pb"' 



10 12 14 16 



A Altarelli, Sdrlinc 
• TokarevAOO 
D TokorevAG<0 



—U- 



P, (GeV/c) 



Figure 6. 

Asymmetry All of 7 + jet production in polarised pp collisions at ^/s 
3 different sets of spin-dependent PDF's (AG >0 ||TT] 



200 GeV for 
as a function of 

transverse momentum p?. The errors indicated are statistical only, based on the expected 
luminosity of RHIC and the properties of the STAR detector. 
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Figure 7. Comparison of the simulated to reconstructed values of the initial-state 
parton kinematics. 
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Figure 8. 

Comparison of the input polarized gluon distribution (line) AG(x) to values extracted 
from the simulated data (points) as a function of the reconstructed x g i uon . The results 
are obtained for spin-dependent parton distribution AG >0 
y/s = 200 GeV. 
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Figure 9. 

Reconstructed values of xAG(x) versus reconstructed values of x g i uon for 3 different 
sets of spin-dependent PDF's (AG >0 [ITTJ, AG <0 III and WW). 
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